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1. Introduction 

This is the write-up of a talk that one of the authors (U.M.) gave on the occasion of the 50th 
Winter Workshop on Nuclear Dynamics in Bormio'. Since this meeting is attended mainly by sci- 
entists working in very different fields, mostly high-energy heavy-ion physics and hadron physics, 
this talk tries to give an impression of what nuclear physics can contribute to neutrino physics, in 
particular with long baseline (LBL) beams. This article thus tries to review some of our work in 
this field for non-experts; as such it draws heavily on our former publications p-p^]. 

Neutrinos are the secondary decay products of pions or kaons produced usually in a high- 
energy proton-nucleus collision. (K2K: 12.9 GeV p on Al and MiniBooNE: 8.9 GeV p on Be) The 
energy-distributions of LBL neutrino beams are thus quite broad (several hundred MeV to several 



GeV). Dedicated experiments, such as HARP [11| and NA61/SHINE [12|, measure the pion and 



kaon production rates and thus provide crucial input for the neutrino beam simulations. For a 
typical nuclear physicist LBL experiments are quite unbelievable: the beam at its source is already 
about 0.5 m wide, at the target it can reach out over several hundred meters. The 'beamlines' are 
long, up to several hundred kilometers, the composition of the beam is known only up to a certain 
extent and, most important, the energy of the beam is distributed in a wide range and is not really 
well determined. Taken all of this together with the fact that typical cross sections for interaction 
with matter are of the order 10^^^ cm^ it is even more remarkable that these experiments have 
actually given us the first look into physics beyond the Standard model. Since neutrino oscillations 
have been verified in different experiments [ pj[ there is now no longer any doubt that neutrinos 
have a non-zero, although quite small, rest mass. At present there are several so-called long- 
baseline (LBL) experiments running with neutrino energies of a few hundred MeV up to several 



GeV [15- 19|. All these experiments aim for a more refined determination of the neutrino mass- 
differences and the oscillation parameters. They also look for fundamental effects such as CP- 
violation in electroweak interactions, the possible existence of so-called sterile neutrinos and cross 



section measurements []20|-|22|]. 

For the determination of the mass-differences the neutrino energy has to be known as is evident 
from the simplified (two flavors only) neutrino oscillation formula 

P(v,^v.) = sin2(20)sin2(^) (1.1) 

which gives the probability that a muon neutrino converts into an electron neutrino Ve over 
the flight path L. Here Am^ is the difference of the squared masses of the two neutrino flavors 
and 6 is the mixing angle. It is seen that with known distance L between source and detector the 
difference of the squared masses is directly proportional to the neutrino energy Ey. The neutrino 
energy is thus a crucial quantity which has to be known on an event-by-event basis. Due to the 
production mechanism, however, only its distribution in the neutrino beam is known. The neutrino 
energy thus has to be reconstructed from the measured energies of the outgoing lepton (for charged- 
current (CC) experiments) and/or the hadrons in the final state of the reaction. For a quasielastic 
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scattering (QE) reaction on a free nucleon at rest the incoming neutrino energy is directly linked 
to the kinematics of the outgoing lepton and is thus known when lepton angle and energy are 
measured. 

Due to practical considerations (difficulties in handling hydrogen targets and to increase the 
count rates) all of these new experiments use nuclear targets (e.g. O, C or Fe). Nucleai" effects 
on the measured inclusive cross sections have been analyzed with the help of Fermi gas models 
that contain the Fermi motion of the nucleons inside the nuclear targets, Pauli-blocking and also 
- roughly - the nuclear binding [23-28]. In addition, nucleons inside nuclei can be collision- 
broadened. For example, Benhar calculates their spectral function within a state-of-the- 
art nuclear many-body theory. Finally, using RPA also corrections of the nucleonic vertex have 



been calculated [B2-37]. Other approaches go beyond the Fermi gas model and in particular a 



plane-wave approximation and calculate QE scattering using more realistic wave functions [|27 



28, 38-43]. Also the so-called scaling analysis of inclusive electron scattering has been used 



for a prediction of neutrino scattering [44-^]. Common to all of these theoretical attempts is 
that they concentrate on QE alone. At energies of a few hundred MeV (as in the MiniBooNE 
and K2K experiments) QE indeed dominates the reaction process. However, even then there are 
higher-energy tails in the energy distributions of the neutrinos so that pion-production events are 
always entangled with QE. All the calculations mentioned are thus not directly comparable with 
'raw' experimental data which always contain some admixture of pion production events. Free QE 
kinematics are then not directly applicable. 

Semi-inclusive processes, in which a nucleon is knocked out from the nucleus, seem to offer 
direct access to QE events. For these processes one needs obviously a very good description of FSI. 
Often, these FSI ai^e completely neglected [ ^ , 41 ] or they are being treated only as absorptive 
interactions, in the framework of an optical model or the Glauber approximation ||2^, 28, |3^, 
42, ^]. In doing so one neglects that collisions of the initially knocked-on nucleons with other 
nucleons in the nuclear target can lead to energy loss, change of direction, charge transfer or to 



multiple nucleon knock-out. These processes can only be described by Monte Carlo models [ 48 



. Contributions to nucleon knock-out come also from resonance excitations. For example, a A 
resonance can first be excited which then decays by a collision with another nucleon thus leading 
to two outgoing nucleons which resemble a QE knock-out event. These 'fake QE' events are very 
important for the experiments since they have a big influence on the energy reconstruction as will be 
demonstrated later on in this article. A quantitative understanding of these events is only possible 
within coupled channel models that contain all the various interaction mechanisms. 

Pion production can not only influence the energy reconstruction, but is essential also for the 
correct interpretation of neutrino oscillation experiments. In particular- neutral pions play an essen- 
tial role in the flavor identification which is essential since every neutrino beam always contains 
admixtures of less dominant flavors. For example, a -produced , which may be produced by 
FSI, can decay to 2y and thus simulate a Ve event in the detector. Up to neutrino energies of about 
1.5 GeV pion production proceeds mainly through the excitation and subsequent decay of the A 
resonance. With increasing energy also higher resonances and non-resonant background can con- 



tribute [ pO[ , pl| , ^]. Singh et al. [ |52| ] have shown that beside medium effects such as Fermi-motion, 
Pauli-blocking and the binding of the nucleons in their potential also more complex in-medium 
effects such as the coUisional broadening of nucleons may influence the width of the A and through 
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it the absorption cross section. Before we started our own work in this field, semi-inclusive observ- 
ables such as pion spectra had been treated only within the old ANP model [^]. 

At high neutrino energies finally deep inelastic scattering (DIS) dominates. Even though nu- 
clear effects have been discussed in this energy region for the NuTeV Experiment [54| (sin^^iy- 
anomaly) most analyses assume that here the nucleus just acts as an ensemble of free nucleons [p^]. 
However, just as there is an EMC effect for electron- and mu-induced reactions on nuclei a similai^ 
effect should be present for neutrinos as well. This is still a field of ongoing research; very recent 
analyzes do not seem to be compatible with the electron-data [56, 57]. 

For a reliable inteipretation of the ongoing neutrino experiments one needs a model that de- 
scribes both the initial neutrino-nucleus interaction and the final state interaction of the produced 
paiticles. Since, because of the always present broad energy distributions in neutrino beams, any 
experiment inherently averages over different reaction types, the model has to be reliable and well 
tested in all three mentioned energy regimes (QE, pion production, DIS). Since none of the works 
mentioned earlier could achieve this goal, special event generators were - quite independently 
from other nucleai^ physics events generators - developed for the simulation of the various neu- 
trino experiments (e.g. NUANCE [||], NEUGEN @, NEUT and GENIE [§l]]). Although 
these different generators are quite different from each other, common to all of them is that they 
all rely on Monte Carlo simulations of the reaction process. These codes contain contributions 
from quasielastic scattering as well as from inelastic processes; for the resonance excitations the 
model of Rein and Sehgal [62] is widely used. They all use a Fermi-gas model for the nucleons 
in the target nucleus, usually with constant Fermi-momentum and constant binding energy. Final 
state interactions are then described by an intranuclear cascade, where the scattering and absorp- 
tion cross sections are adjusted to the special targets used. Very few predictions of these models 
have been compared to data for other reaction types. This is particularly disturbing since there 
exist extensive data for photonuclear and electronuclear experiments on nuclear targets in the rele- 
vant energy regime. Only recently, tests of selected aspects of FSI, i.e. pion absorption, have been 
started [|63[]. Furthermore, the nuclear physics input of these model is often quite doubtful. For 
example, the Rein-Sehgal formfactors fail when applied to electroproduction experiments [ 54 ] and 
in the analysis of the MiniBooNE data even the strength of the Pauli-principle was adjusted to QE 
data [lei. 



2. The GiBUU model 

Over the last 25 yeais we have gained experience in the development of theories, numerical 
methods and codes for the simulation of complete events in nuclear reactions. These methods all are 
based on the Boltzmann-Uehling-Uhlenbeck (BUU) equation, which differs from the mentioned 
neutrino event generators by a consistent treatment of self-energies of all paiticipating particles. 
Fermi-motion in the target nucleus is described by a local Thomas-Fermi approximation and is thus 
closer to empirical distributions than this is normally the case in event generators; the occupation 
in coordinate space follows empirical density distributions. A recent review summarizes some 



of our work with the GiBUU implementation of this method [ \9 ] and gives many more details 



GiBUU essentially factorizes a reaction on a nuclear target into a first, initial interaction (sufficient 
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for inclusive cross sections) and a detailed treatment of FSI (needed for a reliable description of 
semi-inclusive processes). The latter is the particular strength of GiBUU. 

With this method we have initially described the production of quite different particles (pho- 
tons, dileptons, pions, kaons, vector mesons etc.) in heavy ion reactions [ |6^ , 57 1. Later on, we 
have generalized these model studies to more elementary projectiles and have investigated proton- 
induced 



8|], pion-induced []69[ [70||, photon-induced [[71j-[75|] and electron-induced [ [76[ [77| ] par- 
ticle production on nuclei over a wide energy range, using the same theoretical methods and code 
for all these different reactions. In particular, the electron-induced reactions are directly relevant 
for neutrino-induced reactions as they provide a necessary testing ground for the neutrino event 
generators since the initial kinematics is similar and the primary reaction amplitudes are closely 
related. In addition, the hadronic FSI are obviously the same. The electron-induced reactions with 
their large data base are thus ideally suited for a test of the treatment of FSI. 

Indeed, our studies of photon-induced reactions [78, 79] have shown that a detailed, quanti- 
tatively reliable description of the experimental data is possible only if the FSI of the produced 
hadrons is described realistically. In particular, the effects of coupled channels must be taken into 
account which can lead to sidefeeding of the channel under study. In other words, the particle that 
the detector finally - after its traverse through the nucleus - sees must not necessarily be the same 
as the one originally produced in the first interaction. We have shown in studies of electroproduc- 
tion [76] that the importance of such secondary production processes increases with the squared 
momentum transfer Q^. We have also performed extended studies of pion mean free paths and 
their absorption in nuclear targets [79]. Very sensitive to details of the nN interaction are charge 
exchange processes. Their reliable description is thus an important test of the reliability of our 
model for the FSI [0. 

The correct description of the imaginary part of the self-energies of the produced particles, i.e. 
their absorption through interactions with the nucleus, is very important in all these calculations. 
Particles can get 'collision-broadened' in the nuclear medium, i.e. the imaginary parts of their 
self-energies can become large [|^, 81]. It is then a challenge to describe the transport of these 
particles so that they attain their free self-energies when they reach the nuclear surface. This is 
one of the major developments of transport theories in the last decade. We are now in a position 
to transport collision-broadened hadrons consistently and theoretically well founded (so-called off- 



shell transport) [73, B3, 83, 49]. In this point we thus go beyond the simpler approaches in which 



e.g. nucleon resonances are not transported at all, but instead decay immediately; this corresponds 
to a local approximation for the resonance propagators. This point is essential in particular in light 
and medium-heavy nuclei where the change of the density-profile of the target nucleus over the 
propagation-length of the resonance plays an important role. 

At invariant masses higher than about 2 GeV the discrete resonance structure disappears and 
here the region of quark-hadron duality sets in. We, therefore, employ in this regime well-tested 
methods of pQCD in which the total inclusive cross section is described with the help of parton 
distribution functions and the non-perturbative formation of hadrons is encoded in fragmentation 
functions. The latter we obtain from the JETSET (Lund) model which is implemented in the 
code PYTHIA [p4|]. Also this high-energy part we have tested extensively. In particular; we have 
investigated the hadron formation and attenuation in lepton-induced reactions as measured in the 
HERMES and EMC experiments [p5l]. Data from both experiments range from 12 GeV (HERMES) 
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up to 280 GeV (EMC) and can be described simultaneously by our model. These descriptions also 
offer access to the so-called formation time of hadrons during which FSI are suppressed; for details 
see our discussions in [85]. 

Based on all these developments we have extended the GiBUU model over the last few years 
also to neutrino-induced reactions [|l|, ^. The GiBUU model differs from the standard neutrino 
event generators not only by its theoretical foundation. It is, furthermore, the only available method 
that has been tested in comparisons with broad data bases from quite different experiments [^]. 
Nuclear many-body physics is complex and a realistic treatment of FSI requires con^espondingly a 
complex model. This is reflected in the quite sizeable GiBUU code which is available for download 
from [p6|]. There also a description of the numerical implementation can be found; further details 
of the GiBUU model are contained in the review [ffl]. 



3. Results 



As we have discussed already in the introduction one of the essential ingredients for any ex- 
traction of the neutrino masses from oscillation experiments is the neutrino energy. This energy is 
not known a priori in present-days experiments, because neutrino beams are quite broad in energy 
due to their production mechanisms. While at higher energies calorimetric methods may play a 
role, at lower energies (a few hundred MeV to a few GeV) quasielastic (QE) scattering has been 
used to determine the incoming neutrino energy on an event-by-event basis. This method relies 
on an identification of the reaction mechanism (interaction of the neutrino with a single nucleon). 
It also relies on the use of quasifree kinematics that describes neutrino scattering on a single free 
nucleon at rest, thus neglecting any Fermi-motion effects; binding is taken into account only by a 
constant removal energy. 

3.1 The QE-Pion entanglement 

The experimental identification of a special reaction type, in this case QE, is, however, any- 
thing but trivial in an experiment with a broad energy distribution in the incoming beam. There is 
necessarily a mixture of different reaction mechanisms. Up to energies of about 1.5 GeV these are 
mainly QE and pion production. The experiments have thus developed strategies to isolate the QE 
process. In Cerenkov detectors this is done by requiring only the outgoing lepton and the absence 
of any pions. In so-called tracking detectors 1 muon, 1 proton and pions are required for a QE 
event. We call these events 'QE-like'. 

Fig. |l| shows the cross sections for scattering on '^C as a function of neutrino energy both for 
the Cherenkov (/i, Ott) and the tracking detector (^u, 1/7,071) identification methods ^ for QE-like 
events. The solid curve in both cases shows the true CCQE cross section whereas the dashed line 
gives the QE-like cross section. It is seen that the Cherenkov detector identification method leads to 
a cross section that is about 20% higher than the true value; the surplus is due to primary resonance 
or pion production and subsequent FSI leading to pionless final states. For determination of the 
axial mass or for the reconstruction of the incoming neutrino energy this artificial surplus has to be 
removed with the help of an event generator. The opposite is the case for the tracking detector. Here 
the QE-like cross section is about 20% lower than the true value; this is because, first, more than 
one proton may be knocked-out of the nucleus and these events are then not counted as QE-lrke 
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Figure 1: (Color online) Total QE cross section on '^C (solid lines) compared to different methods on how 
to identify CCQE-like events in experiments (dashed lines). The left panel shows the method commonly 
applied in Cherenkov detectors; the right panel shows the tracking-detector method as described in the text. 
The contributions to the CCQE-like events are also classified [CCQE-like from initial QE (dash-dotted) and 
from initial A (dotted lines)]. Experimental detection thresholds are not taken into account ( from |||]). 



and, second, secondary neutrons are not detectable. However, the QE-like sample is very 'clean' 
in that it contains nearly only original QE events. 

A further complication arises because detectors are not perfect but have acceptance thresholds. 
Their effect has been discussed in detail in [^. There it was shown that for realistic detector 
thresholds only about 50 - 70% of all events are observed and that this deficit also affects the 
energy reconstruction. 

3.2 The QE puzzle 

In theoretical calculations the QE cross section is determined by an interplay of vector and 
axial couplings with their corresponding form factors. The vector couplings can be rather well 
determined from electron scattering experiments on the nucleon that work at a fixed energy and 
permit to determine the relevant kinematic parameters, i.e. energy- and momentum-transfer, in each 
event. The corresponding form factors have been shown to have a complicated, non-dipole form 
[p7[]. For the axial couplings the situation is less well determined. Here the data come from electro- 
pion production and older neutrino data on the nucleon or deuterium with large uncertainties. They 
have been analyzed by making a dipole ansatz and then extracting the axial mass from a fit to data. 



The world average for the axial mass parameter is found to be Ma = 1 .026 GeV [ ]88| ] 



It came, therefore, as a surprise when the Mini Booster Neutrino Experiment (MiniBooNE) 



at Fermilab [gOp, which uses a big container filled with oil as a Cherenkov counter, pubhshed its 
results on QE scattering. The analyses of both charged current (CC) and neutral current (NC) high- 
statistics QE events showed a clear excess of cross section over that expected for QE scattering 
from a Fermi-gas model [|^, ^]; see Fig. ^. A similar result had been obtained by the K2K 
experiment |93] that worked with a neutrino flux peaked at a slightly higher neutrino energy (sa 
1 GeV) than the MiniBooNE experiment (ss 0.7 GeV). In both cases the flux distributions are 
rather broad and have a considerable overlap. On the contrary, the NOMAD experiment working at 
significantly higher energies (between about 5 and 100 GeV) observed no such excess of measured 



over expected quasielastic cross section [92|. Both the MiniBooNE and the K2K experiments could 
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Figure 2: (Color online) Flux-unfolded MiniBooNE CCQE cross section per neutron as a function of recon- 
structed neutrino energy. The cross sections are shown along with results from the LSND [ ^ ] and NOMAD 
[ p2[ ] experiments. Also shown are predictions from the NUANCE simulation for a relativistic Fermi gas 
model with two different parameter variations and for scattering from free nucleons with the world average 
Ma value ((f is a parameter influencing the strength of the Pauli -principle). Figure from [|65|]. 



obtain good fits to their data^ in the Fermi-gas model only when the axial mass was considerably 
increased to Ma = 1-23 GeV [ p9| ] or even Ma = 1.35 GeV The sizeable increase in the axial 
mass needed to describe the data cannot be ascribed to deficiencies in the Fermi gas model alone. 
Indeed, Benhar et al. showed that a model based on state-of-the-art nucleon spectral functions 
required an even larger axial mass {Ma = 1.6 GeV) for a fit of the differential data ||3T[]. 



3.3 The pion puzzle 



Recently the MiniBooNE and K2K collaborations have published data on charged [ 94 1 and 



neutral [95] pion production in CC neutrino scattering. Generally, these experiments report cross 
sections which are noticeably higher than those expected from any conventional theoretical ap- 
proach @ |]. 

As an illustration we present our calculations on a CH2 target and compare the results with the 



data from the MiniBooNE experiment [ ]94[ , |95| ]. Integrated cross sections versus neutrino energy for 
the charged current l7r+ and Itt" production are shown in Fig. ^ As in the MiniBooNE experiment, 
the one pion events are defined as "observable one pion production", i.e. events with one pion of 
a given charge and no other pions in the final state, regardless of which particles were produced in 
the initial neutrino vertex. 

The left upper panel in Fig. ^ shows the results for 71+ production. Comparison of the curves 
with and without FSI shows that the FSI do not change the energy-dependence of the curves. The 
curve with FSI (solid line) lies 20% below the curve without FSI (dash-dotted line). This reduction 
is mainly caused by the A absorption through — )■ A^A'^ scattering. Charge exchange n^n n^p 
also depletes the dominant l7r+ channel. The side feeding from the reverse process gives minor 
relative contribution because the initial tt" production cross section is around 5 times lower. For 
the Itt" cross section in the left lower panel the same side feeding processes increase the cross 
section. This is, however, compensated by absorption, charge exchange to the 71^ channel through 



^It must be noted that neither the cross section nor the energy in Fig. ghave been directly measured; both involve 
some model dependence. 
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Figure 3: (Color online) Integrated cross section for MiniBooNE l;r+(upper panel) and l7r''(lower panel) 



CC production versus reconstructed neutrino energy. Data are from ||94[ The panels on the right show 
the composition of the calculated cross sections. Figure from [P] 



K n K p and other channels such as kN — ?• LK,AK. The overall effect is nearly the same cross 
section with and without FSI. 

The right upper panel shows the origin of the l7r+ events. Most of them (dash-dotted line) 
come from initial A resonance production and its following decay. Some events (dotted line) are 
background ones. A minor amount comes from the initial QE vertex (long-dashed line), which is 
only possible due to FSI, when the outgoing proton is rescattered. Here the main contribution is 
from the pN — > A'^'A — > N'n" K reaction. The right lower panel shows the origin of the iTl" events. 
Here the background processes and the FSI play an even bigger role. 

We now turn to a discussion of a comparison with experiment. Already the curves without 
any FSI lie considerably {k, 25%) below the data, those with FSI included (solid line) are nearly 
by a factor of about 1.5 (at 1 GeV) below the experimental data; at higher neutrino energies these 
discrepancies become even larger since the experimental cross sections rise more steeply with 
energy than the calculated pion cross sections. At 2 GeV the discrepancy between the results with 
FSI and the data amounts to a factor of « 2. For 1 7r° production one observes the similar result: 
the data are considerably higher than our calculated values. 

Fig. ^ shows our calculations for the kinetic energy distribution of the outgoing pions. As for 
the other distributions, our calculations with FSI are lower than the experimental data by a factor 
of 1.6-2. 

A distinctive feature of our calculations is the result that FSI significantly change the pion 
spectra. The lowering of the cross section for Tj^ > 0.12 GeV is a direct consequence of the pion 
absorption through TlN — >■ A with the following AA'^ — )• A^A'^. Pion elastic scattering in the FSI also 
decreases the pion energy, thus depleting spectra at high energies and accumulating strength at low 
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Figure 4: Distribution of (left panel) the outgoing 7r+ in their kinetic energy; (right panel) the outgoing K in 
their absolute value of the 3-momentum for the MiniBooNE CC neutrino scattering. Data are from 
figure from [^] . 



energies. For production an additional increase of the cross section at lower energies comes 
from the side feeding from the dominant ;r+ channel as discussed above. The change of the shape 
of the spectra is similar to that calculated for neutral current l7r° production in [||]. 

The predicted shape of the pion spectra is due to well-known ;r — A'^ — A dynamics in the 
medium is the same as that observed experimentally in (7,71") production on nuclear targets [78 
97]. Its absence in the neutrino data is, therefore, hard to understand |98| ]. 



A source of uncertainties in these comparisons is the determination of the neutrino energy 
which is done using quasi-free kinematics for on-shell Delta production on a nucleon at rest. One 
more point is that the experimental data are presented after being corrected for the finite detection 
thresholds for muons, pions and nucleons. This may introduce an additional dependence on a 
particular neutrino event generator used in a given experiment 

3.4 Energy Reconstruction 

We now discuss the quality of the energy reconstruction which is based on applying quasifree 
kinematics of true QE scattering to QE-like events and neglecting Fermi motion Fig. ^ shows 
the distribution of reconstructed energies, obtained in a GiBUU simulation, for a fixed incoming 
neutrino energy of 1 GeV. The distribution is clearly affected by the two effects discussed above: 
Fermi motion leads to a broadening of the neutrino energy around the incoming energy; this ac- 
counts for the broadened peak at 1 GeV. Its width of about 16% is determined by Fermi motion 
alone and is thus always present; this defines the lower limit for any energy reconstruction via QE 
scattering. In addition the reconstructed energy distribution exhibits a clear bump at lower energies; 
this is caused by an initial pion production. The bump does depend on the pion detection threshold 
and amounts to about 15% of the true QE peak height at a realistic detection threshold of about 100 
MeV pion kinetic energy. Taking this lower-energy bump into account raises the rms energy-width 
to about 22 % (for a more detailed discussion see [||]). A qualitatively similar result has also been 



obtained by Tanaka (see Fig. 3 in [100]) using the GENIE event generator; there the lower-energy 
bump is even more pronounced. 

The T2K experiment has recently reported the first experimental observation of electron neu- 
trino appearance from a muon neutrino beam [p9|]. At the flux maximum the energy of is 600 
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Figure 5: (Color online) Distribution of the reconstructed neutrino energy using quasifree kinematics on a 
nucleon at rest for QE and for static A excitation for fij^^' = 1 GeV. Shown is the reconstruction based on 
the CCQE-like sample (before and after fsi and Cherenkov assumptions) and based on the CCl7r+ sample 
(before and after fsi). Notice that the dash-dotted and dotted curves partially overlap (from |M|). 
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Figure 6: (Color online) NC induced single- ;r'^ production cross sections on ^^O averaged over the T2K 



flux and multiplied by the misidentification probability (dashed line, taken from Ref. [ 101 ]) as a function of 
the pion momentum (from [^). 



MeV. At this energy the error in the reconstructed energy, according to Table I in ||8]|, amounts to 
about 21% for the Cherenkov detector. Fig. |I| shows that at 600 IVIeV the QE events are fairly clean, 
but that the admixture of pion production events rises significantly, already up to 1 GeV, reflecting 
the pion production threshold. 

A crucial problem in any such experiment with Cherenkov counters is that the decay photons 
of neutral pions can be misidentified as electrons. We therefore, investigate in Fig. ^ the probability 
that a misidentified %^ is counted as a appearance event. The probability that a %^ cannot be 



distinguished from is given by the dashed line (taken from Fig. 2 of Ref. 1 101 1). The solid line 
shows the weighted cross section averaged over the T2K flux and calculated using the elementary 
ANL data as input. The total cross section for misidentified events is now 0.09 • 10^^^ cm^ and 
thus 26% of the true pion events. 



3.5 Many-body interactions 

All the results discussed so far rely on the impulse approximation. It has, however, recently 
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been pointed out that the identification method of the MiniBooNE experiment allows for a signif- 



icant amount of non-QE many -body excitations in the QE-Iike cross section IIO4 p7| , |103| ]. 
These excitations are connected with more than one outgoing nucleon in the initial, primary inter- 
action and can be caused by the interaction of the incoming neutrino with more than one nucleon. 
That such reactions contribute to the inclusive electron scattering cross sections on nuclei is quite 
well known since about 20 years when it was shown that such events contribute significantly in the 
so-called dip region between the QE peak and the A resonance. For a more extended discussion of 
the literature see [10|. 

These so-called 2 particle - 2 hole (2p-2h) effects can indeed account for the disagreement 
between the data in Fig. ^ and the calculations using the world average value of about 1 GeV for 
the axial mass and are not incorporated in the energy reconstruction discussed in the last section. 
Since for the 2p-2h excitations the quasifree kinematics formulas used by experiment for the energy 
reconstruction do not apply, the actual energy uncertainty may even be larger than discussed before. 
There are two features of of these 2p-2h contributions that affect the energy reconstruction. First, 
at forward angles, where the cross section is largest, the 2p-2h contributions are largest at small 
T^, below the QE peak. When analyzing such events with the help of the one-body expression 
this leads to a lower reconstructed neutrino energy than the true one. Second, while QE scattering 
is strongly forward peaked, the 2p-2h events are fairly fiat (within a factor of 2) in lepton angle 
(see [|To|]). The relatively strong yield at backward angles will lead to a larger neutrino energy, in 
particular for intermediate muon energies. Since both effects are present we expect a fairly flat 
behavior of the 2p-2h contribution to the energy reconstruction. 




Figure 7: (color online) Cross section per energy per nucleon for QE-like processes as defined by the 
MiniBooNE experiment (0 pions and any number of nucleons in the final states) versus the reconstructed 
neutrino energy for the true energy of the incoming neutrino of 0.5, 0.7, 1 .0 and 1 .5 GeV. Figure from [[lO|]. 

We have recently implemented the 2p-2h excitations into GiBUU |1C]. In Fig. ^ we plot 
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the distribution of the reconstructed neutrino energy obtained using the MiniBooNE reconstruc- 
tion method with Eb = MeV (which is a typical binding energy in the GiBUU code, as op- 
posed to the value —34 MeV used by the MiniBooNE) for the fixed true neutrino energies Sj,™ = 
0.5, 0.7, 1.0 and 1.5 GeV. The main contribution to all QE-like events is given by the true QE 
events, which shows a prominent peak around the real energy. The peak is approximately symmet- 
ric and has a width of about 0.1 GeV: this broadening is caused by Fermi motion. For A— induced 
events (as already discussed in ^ and in the previous section) the distribution is not symmetric, 
with a broad peak at lower energies. For a more detailed discussion of these effects see [^, ^. 

For all energies the broad 2p-2h contribution is peaked below the true energy. For the lower 
neutrino energies of 0.5 and 0.7 GeV most of the total distortion is caused by 2p-2h processes 
because at this low energy there is only little A excitation. For the true neutrino energies of 1 .0 and 
1.5 GeV, on the other hand, the contributions of 2p-2h events are comparably small because here A 
excitation plays a major role and because their strength is distributed over the whole energy range 
from 0.2 to 1.8 GeV, with a flat maximum around 1/2 of the true energy. Both the 2p-2h effects as 
well as the A excitation lead to a shift of the reconstructed energy towards smaller values, or, vice 
versa, for a given reconstructed energy the true energy always lies higher than the reconstructed 
one. The effect is most pronounced at lower true energies and becomes smaller at higher energies. 
Both of these results agree with a recent analysis by Martini et al. 1 104 1. 




Figure 8: (color online) Event distribution (i.e. the cross section multiplied by the neutrino flux) for QE- 
like processes as defined by the MiniBooNE experiment versus the real (left picture) and reconstructed 
(right picture) neutrino energy. The thin solid and dot-dashed curves repeat the distributions for true and 
reconstructed energy in both figures. Figure from [p^. 

This is also seen in Fig. ^ which shows the event distribution (that is the cross section mul- 
tiplied by the neutrino flux) for QE-like processes that would be observed by the MiniBooNE 
experiment. All events with zero pions in the final state are identified as QE-like ones. Since in the 
numerical simulation we know the true neutrino energy and the original interaction vertex, we can 
predict the origin of each event and its energy distribution for a given true neutrino energy. Only 
parts of these events are genuine QE events. The same QE-like events can originate from initial 
A production, higher resonance production (not shown in the figure), pion-background events (not 
shown in the figure) or 2p-2h processes. The event distribution versus the true energy is shown by 
the curve labeled 'true Ey . Comparing it to the one versus the reconstructed energy (rec Ey) one 
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sees that there is a systematic distortion in energy reconstruction: the peak of the event distribution 
is shifted now by about 100 MeV to lower energies with this shift becoming smaller for the larger 
energies above the peak. Qualitatively this same behavior was already observed and discussed for 
the K2K and the MiniBooNE flux in [||] as a consequence of pionic excitations. Here, now also 
the 2p-2h contribution is included in this analysis; it leads to the noticeable downward shift of the 
reconstructed energy curve on the low-energy side of the maximum where the pion production is 
still small. 



4. Summary 

In summary, understanding the reaction mechanism of neutrinos with nuclei is mandatory for 
a determination of neutrino oscillation parameters. Many of the relevant cross sections are experi- 
mentally not determined and have to be provided by theory. Here neutrino physics can only benefit 
from a close contact with photonuclear physics. This is true for the elementary cross sections, but 
it is also true for the all-important final state interactions. The latter are independent of the primary 
interaction process and can - and must - be checked with a multitude of other nuclear reactions. 
Among nuclear physicists working in heavy-ion physics and electron-, photon- or hadron-induced 
reactions on nuclei, there exists lots of experience and knowledge that still has to find its way into 



the neutrino community [105|. In particular the all-important neutrino event-generators can only 
benefit from checks with the help of relevant photonuclear data which are readily available. 

One important example for the influence of nuclear effects can be found in the determination 
of the neutrino energy which has to be known for the extraction of the neutrino mass- and mixing- 
parameters. Using quasielastic scattering as the main method to determine this energy requires a 
good identification of this reaction process. This involves the separation of 'clean' QE processes 
from resonance-excitations and many-particle interactions. This is illustrated very nicely in the 
results of the MiniBooNE experiment that have found widespread interest in the community over 
the last few years. The large excess of QE-like events observed by MiniBooNE over those cal- 
culated by commonly used event generators with the world-average value for the axial mass of 
around 1.0 GeV has led to a multitude of partly contradictory models that all aim to describe the 
same data starting from quite different initial Ip-lh or 2p-2h interaction mechanisms. Since broad 
energy-band neutrino experiments necessarily contain an averaging over quite different reaction 
mechanisms the inclusive data alone do not allow for an unambiguous experimental verification of 
any particular model. We have discussed that the energy reconstruction is affected by the presence 
of 2p-2h interactions mainly at the lower neutrino energies while the A excitations become more 
essential at the higher energies. 2p-2h interactions lead - as the A excitations do - to a downward 
shift of the reconstructed energy. Their contribution to the reconstructed energy distribution is 
fairly flat over a wide reconstructed energy range for the higher true energies. 

In this talk we have given a short discussion of the relevant ingredients for an understanding 
of neutrino-nucleus data. Quasielastic scattering, pion production and many-particle interactions 
all have to be well controlled. The presently available inclusive data all rely for their understanding 
on event generators. A closer check can, however, probably come only from more exclusive data, 
like those on pion production and nucleon knock-out spectra. Nuclear physics can contribute a lot 
to understanding these data. 
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